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' Abstract. We present an introduction to coined quantum walks on regular graphs, 

^SJ , which have been developed in the past few years as an alternative to quantum Fourier 

I ' transforms for underpinning algorithms for quantum computation. We then describe 

I our results on the effects of decoherence on these quantum walks on a line, cycle and 

hypercube. We find high sensitivity to decoherence, increasing with the number of steps 
in the walk, as the particle is becoming more delocalised with each step. However, the 
\ effect of a small amount of decoherence can be to enhance the properties of the quantum 

walk that are desirable for the development of quantum algorithms, such as fast mixing 
times to uniform distributions. 
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Quantum walks are based on a generalisation of classical random walks, which 
, have found many applications in the field of computing. Examples of the power 

' of classical random walks to solve hard problems include algorithms for solving 

, fc-SAT [1], estimating the volume of a convex body [2], and approximation of 

^ • the permanent [3]. They are a subset of a wider model of computation, cellu- 

lar automata, which have been proved universal for classical computation. The 
P5 ' utility of classical walks suggests that extending the formalism to the quantum 

, regime may assist the new field of quantum information processing in generating 

further quantum algorithms. Similarly to the classical case, it is also possible to 
define the notion of quantum cellular automata, whose equivalence to quantum 
Turing machines has been shown [4]. 

Most known quantum algorithms are based on the quantum Fourier trans- 
^ • form, for an introduction to quantum computing and algorithms see, e. g., [5, 6]. 

Quantum versions of random walks provide a distinctly different paradigm in 
which to develop quantum algorithms. Very recently, two such algorithms have 
been presented. Shenvi et al. [7] proved that a quantum walk can perform the 
same task as Grover's search algorithm [8], with the same quadratic speed up. 
Childs et al. [9] describe a quantum algorithm for transversing a particular graph 
exponentially faster than can be done classically. This exponential speed up is 
very promising, though the problem presented is somewhat contrived. 

In fact, several possible extensions of classical random walks to the quantum 
regime have been proposed [10-12], however, here we will only treat the discrete 
time, coined quantum walks [13], subsequently these are referred to simply as 
"quantum walks" . Before introducing quantum walks, it is helpful to review the 
properties of classical random walks. This is followed by an overview of quantum 
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walks on a line, A''-cycle, and hypercube. Section 2 presents our results on the 
effects of decoherence in these quantum walks. 

1.1 Classical Random Walks on Graphs 

The discrete space on which a random walk takes place can most generally be 
described as a graph G{V, E) with two components, a set of vertices V, and a 
set of edges E. An edge may bo specified by the pair of vertices that it connects, 
e = (vj, Vo). The graph is undirected when (u,, Vq) G E \S.{vo, Vi) G E. The second 
essential feature of a classical walk is the (time-independent) transition matrix 
M, whose elements Mij provide the probability for a transition from vertex Vi to 
vertex Vj . These probabilities are non-zero only for a pair of vertices connected 
by an edge. 

My ^0 iff e={vi,Vj)eE. (1) 

The walk is unbiased if the non-zero elements of M arc given by My = j-, where 
di is the degree of the vertex Vi (the number of edges connected to Vi). A graph 
is called regular (d-regular) if all vertices have equal degree (rf). The state of a 
classical random walk at a given time t is described by a probability distribution 
P{v,t) over the vertices v &V. This distribution evolves at each time step by 
application of the transition matrix M, 

P(w,t) = M*P(w,0) . (2) 

A number of features of these classical walks are worthy of note for later 
comparisons. If G is connected (every pair of vertices have a path linking them 
via a sequence of edges), then the walk tends to a steady state distribution tt 
which is independent of the initial state P(v, 0). Further, if the graph is regular 
then the limiting distribution is uniform on all the vertices. An exception arises 
in the case of periodic random walks, however, a "resting probability" for the 
walk to remain at the current vertex may be added which breaks the periodicity 
and restores the usual convergence. The rate of convergence to this limiting 
distribution may described in a number of ways, here we choose the mixing 
time, 

Mf = min {T\\/t > T : | \P{v, t) - tt\ |tv < e} • (3) 
The measure used here is the total variational distance, 

IIPI -P2||tv = \pi{Vi) ~ P2{Vi)\ . (4) 

viev 

The mixing time thus gives a measure of the time after which the distribution 
is within a distance e of the limiting distribution and remains at least this close. 
An alternative measure that is useful in classical algorithms, for example 3- SAT 
[1], is the hitting time. This is defined for a pair of vertices vq, vi as the expected 
time at which a walk starting at vq reaches v\ for the first time. 
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1.2 Criteria for a Quantum Walk 

Given a d-regular graph, G{V, E), the associated Hilbert space may be defined 
as 

Hv = span{\v,)}Z[ ■ (5) 

As aheady mentioned, several extensions of classical random walks to the quan- 
tum regime have been proposed, there is no unique way to do this. To guide 
us, wc note there arc several properties of classical random walks and quantum 
systems that we would like such quantum walks to have, if possible. The three 
desirable properties of the classical transition matrix M are, 

• locality: My 7^ iff e = {vi,Vj) G E, i. e., transitions are only between 
vertices directly connected by an edge 

• homogeneity: if e = {vi,Vj) G E, |My| = ^, i. e., equal probability of 
transition to any neighbouring vertex 

• time independence: the current step does not depend in any way on previous 
steps of the random walk (Markovian) 

The quantum transition matrix U, equivalent to M, should have similar prop- 
erties to M. In addition, it should be 

• unitary 

because all pure quantum processes are unitary. As proved by Meyer [14, 15], this 
requirement of unitarity is incompatible with the three prior properties, except 
in very special circumstances that don't produce interesting quantum dynamics. 
In order to generate a non-trivial quantum evolution, one or more of these con- 
straints must be relaxed. The continuous time quantum walks proposed by Farhi 
and Gutmann [10], is non-local, in the sense that there is a small probability of 
the particle moving arbitrarily far away in a given unit time interval. The quan- 
tum cellular automata of Meyer [15, 16] are not homogeneous, in that the full 
dynamics is specified over two time steps rather than one. One may make the 
quantum walk non-unitary by measuring the particle at each step, but this sim- 
ply reproduces the classical random walk. Finally, one may consider relaxing the 
time independence condition slightly, and this is what the coined quantum walks 
effectively do. For completeness, we note that there is an equivalent formulation 
of the coined quantum walk in terms of a simple quantum process on a directed 
graph, derived from the original undirected graph, that is due to Watrous [17]. 

1.3 Coined Quantum Walks 

By analogy with the classical walk, in which one pictures flipping a coin at each 
time-step to determine which edge to leave the current vertex by, interesting 
quantum results may be obtained by the introduction of an explicit (quantum) 
coin. In addition to the particle, with associated Hilbert space Hy as given 
in (5), there is an auxiliary quantum system, the coin, with a Hilbert space 
of dimension d (the degree of the graph) He = C'^- The transition matrix is 
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then defined as a unitary matrix U acting on the tensor product of these two 

spaces. This unitary is constructed from two separate operators, a "coin flip" 
and a conditional translation. For an unbiased walk, the coin flip operator, C is 
a unitary matrix whose elements all have equal modulus in the computational 
basis of the coin. This adds significant new degrees of freedom to the system, as 
the relative phases of these elements may be chosen arbitrarily. The conditional 
translation T moves the particle along an edge to an adjacent vertex determined 
by the state of the coin. The evolution of the system from an initial state |V'(0))cv 
to |'0(i))cv S'ftsr t steps is thus given by, 

IV'(i))cv = [T-(C®I)]*|V'(0))cv (6) 

The coin can also be thought of as forming a d-statc quantum memory from one 
step to the next, allowing a much wider range of dynamics to be fully reversible 
(a necessary property of unitary evolution). 

Unitary evolution is also completely deterministic, so the choice of initial 
condition is never "washed out" , rather, it plays a key role in determining the 
outcome of the quantum walk. Unitarity also means that the joint state of the 
coin and particle can never reach a steady state. Even the induced probability 
distribution over the nodes obtained by tracing out the coin, 

P{V, t) = y{v\Tvc mt)}cvm)\] , (7) 

does not converge to a long-time limit. However, it is possible to define a time- 
averaged probability distribution that does reach a steady state, 

pKT)=i:^. (8) 

Operationally, this is easy to produce, it is the distribution obtained by sampling 
the particle location at some time t uniformly selected at random from < t < T. 
Using this distribution it is possible to define a quantum mixing time similar to 
that in (3), 

= min {T|Vt > T : | \P{v,t) - tt\ |tv < e} • (9) 

The second measure discussed for classical random walks was the hitting 
time. This can also be extended to quantum walks, but as before will require 
some modification. In the classical case it is possible to measure the location 
of the particle at each time step to ascertain whether or not it has reached the 
desired vertex. The equivalent quantum measurement disturbs the system; if 
a complete projection onto the vertex Hilbert space is performed at each step 
then all coherences are lost and a classical distribution results. There are two 
approaches to this problem that have been used in the literature [18]. It would 
be possible to wait until a chosen time, T, and then perform a full measurement 
over all the vertices. If the probability of being at the desired vertex is greater 
than some value p which is bounded below by an inverse polynomial in the size 
of the graph, it is said that the walk has a (T,p) one-shot hitting time. (This 
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lowerr bound ensures that standard amplification procedures can efficiently raise 

the success probability arbitrarily close to unity.) The second alternative is to 
perform a partial measurement at each time step. Projecting onto the subspaces 
given by the desired vertex v, IPy = \v){v\, and its orthogonal complement, 
JPg = I — \v){v\ will halt the walk once the vertex v is reached. A walk has 
a {T,p) concurrent hitting time if IPy is measured with probability greater or 
equal to p, which again must have as a lower bound an inverse polynomial in 
the graph size, in a time T. Both these quantities have been shown to exhibit an 
exponential speed up over the classical case for a walk between opposite corners 
of a hypercube [18]. 

Explicit solution of a quantum walk has proved to be a difficult problem. In 
certain cases standard techniques from classical graph theory have been applied 
with some success, namely solution in the Fourier space of the problem [13, 19], 
and also the technique of generating functions [20] . Fourier solutions are possible 
when the graph is of a particular form, known as a Cayley graph. Any discrete 
group has an associated Cayley graph, in which the elements of the group form 
the vertices and the edges are placed by choosing a complete set of generators, 
gi and placing an edge between the vertices a and 6 iff, a = g^b for one of the 
gi. This produces an undirected graph results if the group is Abelian. We will 
next consider three specific cases of such graphs for which analytical solution 
has proved possible. 



1.4 Coined Quantum Walk on a Line 

Consider the case when the graph consists of an lattice of points at integer 
positions on an infinite line. The Hilbert space of the particle is represented by 
the integers, 

Hv = span{|a;) : x G Z} . (10) 

At each time interval the particle can move one step to either the left or the 
right, the coin thus requires only a two dimensional Hilbert space which can 
conveniently be written as. 

He = {\R),\L)} . (11) 

The conditional translation T thus has the following action on the basis states, 

T\R,x) = \R,x + l) , T\L,x) = \L,x-l) . (12) 

The evolution of the probability distribution induced on the lattice may now be 
compared with that of the analogous classical walk on a line. It is sufficient to 
us the Hadamard operator as the coin flip, 

" = ;7l(l-i). ('^) 

all others are essentially equivalent [19, 20]. We also choose the initial state to 
be 

\m) = ^m+i\L))®\o). (14) 
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Fig. 1. A schematic representation of tiie first four steps of tlie evolution of tlie 1-D 
walk with a Hadamard coin. Circles show the particle positions and arrows show the 
state of the coin, with relative magnitudes indicating the amplitude of the components 
of the wavefunction 



This choice produces a symmetric distribution, as in the classical case. Since H 
and T contain only real elements, there is no interference between the part of the 
walk due to the initial state \R) and that due to i\L). Thus the natural (in terms 
of basis states) outcome of a quantum walk is actually biased, unlike a classical 
random walk, a basic demonstration of the initial condition affecting the entire 
outcome of the walk. Four steps of this evolution are shown schematically in 
Fig. (1) and the respective probability distributions are shown in Fig. (2) after 
100 time steps. The classical walk forms a binomial distribution, with a mean of 
and standard deviation ac = Vt. Figure 2 displays several of the interesting 
features of a quantum walk. The central interval [— 2\/T, 2\/T] is essentially 
uniform in the quantum case, with oscillating peaks outside this region up to 
[—T/y/2, T/y/2]. This is due to interference between the large number of possible 
paths to each point in this range. 

The quantum walk on a line has been solved exactly [19, 21] using both 
real space (path counting) and Fourier space methods. The solutions are com- 
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Fig. 2. A comparison of the probability distribution of a classical and coined quantum 
walk on a line, after 100 time steps. Only even points are plotted as both functions are 
zero on the odd grid points 



plicated, mainly due to the "parity" property, namely, that the solutions must 
have support only on even(odd)-numbered lattice sites at even(odd) times. The 
moments can be calculated for asymptotically large times, for a walk starting at 
the origin, the standard deviation is 

1 

Til ^ 

The standard deviation ctq is thus linear in T, in contrast to ^/T for the classical 
walk. 

Quantum walks on a line have now been studied in considerable detail. Dis- 
cussions of absorbing boundaries have been given [19, 20, 22], with applications 
to halting problems in mind. Extensions to multiple coins have been made by 
Brun et al. [23, 24]. However, though useful for understanding the basic prop- 
erties of quantum walks, the walk on a line is too simple to yield interesting 
quantum problems for significant algorithms. 



lim 

T^oo 




1.5 Coined Quantum Walk on a AT-Cycle 

If periodic boundary conditions are applied, instead of a walk on an infinite line, 
a closed walk on a TV-cycle is obtained. For this system, which wraps round 
on itself, the standard deviation is inappropriate and mixing times must be 
considered. The classical result is well known, 

= o{N^/e) . (16) 

Here the time average of the probability distribution, P(v,T) has been used 
in (3) for ease of comparison with the quantum case. The usual mixing time 
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defined by (3) using P{v,T) scales as = o{N'^log{l/s)). The scaling of e 
is not important for the quantum-classical comparisons wc will do here, more 
details can be found in [13]. The equivalent quantity for a quantum walk was 
bounded above by Aharonov et al. [13], and shows a quadratic speed-up, 

Bounds have also been established in the case of more general graphs and it is 

conjectured that mixing times can be improved at most polynomially by quan- 
tum walks [13]. Numerical studies, see [25], suggest the form of the quantum mix- 
ing time is actually = O {N/e), and explain why tighter analytical bounds 
are hard to obtain. 

1.6 Coined Quantum Walk on a Hypercube 

A hypercube of dimension is a Cayley graph with 2^ vertices labelled by the 
bit-strings of length N. The Hilbert space for this graph is spanned by the basis 
states \x) . X G [0, 2''^]. Each vertex has degree A, so the state space of the coin is 
Tic = C^^. The basis vectors of this space are denoted by |a), a G [0, A]. These 
states correspond to the N vectors \ea) where Co is the vector with all zeroes 
except for a single one in the a**^ position. The translation operator for the walk 
on a hypercube can then be defined as, 

T\a,x) = \a,x®ea) . (18) 

Any A X A unitary matrix with all elements of unit modulus may be used as an 
unbiased coin, however, this is not the most natural choice given the symmetries 
of the hypercube. Instead, a biased coin has been selected, which distinguishes 
the edge along which the particle arrived at the vertex from all the others [26]. 
This "Grover" coin acts thus, 

g|«) = |EI^)-I«)- (19) 

h 

The full evolution for a single time step of this walk is then given by U = 
T- (G0I). 

The quantum walk on the hypercube can be solved analytically by mapping 
it to a walk on a line with a variable coin operator, provided the symmetry of the 
hypercube is maintained throughout. To preserve the symmetry between particle 
states with equal Hamming weights, defined to be the number of I's in their bit- 
string, the initial state must be chosen to be localised on the hypercube at |0) 
with the coin in an equal superposition of all its states, |V'(0)) = X^o 1^) ® jO)- 
With this ansatz it has been shown by Moore and Russell [26] that the walk has 
an exact instantaneous mixing time, when the distribution transiently becomes 
exactly uniform, at T = irN/'l. More importantly, Kempe [18] found the first 
exponential gap between classical and quantum walks for hitting times on the 
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hypercube. Both the concurrent and one-shot hitting times for a walk to travel 

from |0) to |2^) are foimd to be ttN/2, with small error probability. This is 
exponentially faster than the classical hitting time, T = 2^~^. This is not a 
true quantum speed up, since there are more sophisticated classical algorithms 
than a random walk to reach the opposite corner that exploit the symmetry of 
the hypercube. Nonetheless, it is important as the first indication that quantum 
walks have the potential for driving quantum algorithms with exponentially 
faster, later confirmed by Childs et al. [9] on a more random graph. 

2 Decoherence in Quantum Walks 

We will now study each of the three systems just introduced in the presence 
of decoherence. A simple model of decoherence is chosen here, at each time 
step of the quantum walk a measurement is made in the computational basis 
with a probability p. We consider three cases, where the measurement is over 
only the coin degrees of freedom, only over the particle states, or is a complete 
measurement of both. Given a unitary transform U = T • (C (8) I) for each step 
of the walk as described by (6), the effect of this decoherence model can be 
described as a discrete master equation, 

Q{t+l) = {l-p)Ug{t)U^ +pJ2^iUQ{t)U^IPi ■ (20) 

i 

The summation runs over the dimensions of the Hilbert space on which the 

decoherence occurs, either the coin He, the particle Hv or both He ^ Hy- The 
projectors JPj are defined to act in the computational basis. When p = the 
ideal quantum walk is obtained, and for p = 1 when a measurement is made at 
each time step it produces the classical random walk. 

The numerical work summarised in the remainder of this paper has been 
presented more fully in [25]. 

2.1 Decoherence in a Quantum Walk on a Line 

Since there is an experimental proposal to implement a quantum walk on a line 
in an optical lattice [27], as well as the three examples for JPj given above, we 
considered the likely form of experimental errors, and also modeled the effect of 
an imperfect Hadamard on the coin. The Hadamard operation may be considered 
to be a "rotation" about the computational basis by 7r/4, (actually a rotation 
and reflection since det(Rot(^)) = — 1), 

Rot(^)=f^i"f) ' (21) 

The error model used in this case consisted of a Gaussian spread of standard 
deviation ^/pn/'i about the ideal value of ^ = 7r/4. 

All types of decoherence model produce the same general form for the decay 
of (Tp{T) from the quantum to the classical value, with small differences in the 
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Fig. 3. Standard deviation ap{T) of the particle on a line for different models of deco- 
herence, for T = 100 time steps 



rates, as shown in Fig. 3. The slope of crp{T) is finite as p ^ and zero at 
p = 1. We have calculated (Jp{T) analytically for pT <C 1 and T 1 for the case 
where JPj is the projector onto the preferred basis {\a, a;)}, i. e., the decoherence 
affecting both particle and coin. Details are given in [28]), the result is 

ap(T) < a{T) 

The first order dependence is thus proportional to pT , so the sensitivity to 
decoherence grows linearly in T for a given decoherence rate p. 

There are interesting differences in the shape of the distribution of the particle 
position for each of the types of decoherence. The decoherence rate that gives 
the closest to uniform distribution has been selected and plotted in Fig. 4, along 
with the pure quantum and classical distributions for comparison. When the 
particle position is subject to decoherence that tends to localise the particle in 
the standard basis, this produces a highly uniform distribution between ±T/ \/2 
for a particular choice of p. The optimal decoherence rate p„ can be obtained 
by minimising the total variational distance between the actual and uniform 
distributions, 

i/(p,T)= ||P(x,p,T)-P„(T)||tv = ^|P(x,p,T)-P„(T)| , (23) 

X 

where P{x,p,T) is the probability of finding the particle at position x after 
T time steps, regardless of coin state [compare (7)], and Pu{T) = V^/T for 
-T/V2 < X < T/V2 and zero otherwise. The optimum decoherence rate de- 
pends on the number of steps in the walk, it can be determined numerically that 
PuT ~ 2.6 for decoherence on both the particle and the coin, and PuT ~ 5 for 



i;^ + ^(i-i/v^) + o(pM/r) 



(22) 
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particle position 

Fig. 4. Distribution of tiie particle position for a quantum walk on a line after T = 
200 time steps. Pure quantum (dotted), fully classical(dashed), and decoherence at 
rate shown on part of system indicated by key (solid). Uniform distribution between 
—T/y/2 < X < T/\/2 (crosses) also shown 

decoherence on the particle only. These differences in the quality of the uniform 
distribution are independent of p and T, and provide an order of magnitude (0.6 
down to 0.06) improvement in v over the pure quantum value. Decoherence just 
on the coin does not enhance the uniformity of the distribution, as Fig. 4 shows, 
there is a cusp at a; = 0. 

2.2 Decoherence in a Quantum Walk on a AT-Cycle 

We now consider a walk on a TV-cycle subjected to decoherence. There is an 
experimental proposal for implementation of a quantum walk on a cycle in the 
phase of a cavity field [29], in which further aspects of decoherence in such 
quantum walks are considered. Recall from Sect. 1.5 that the pure quantum 
walk on a cycle with N odd, is known [13] to mix in time < O(A^logA^) if a 
Hadamard coin is used. The quantum walk on a cycle with N even does not mix 
to the uniform distribution with a Hadamard coin, but can be made to do so by 
appropriate choice of coin flip operator [30] . Under the action of a small amount 
of decoherence, the mixing time becomes shorter for all cases, typical results are 
shown in Fig. 5. In particular, decoherence causes the even-TV cycle to mix to 
the uniform distribution even when a Hadamard coin is used. The asymptotes in 
Fig. 5 for N even and decoherence on the coin only, for p < 2/N, are well fitted 
by epM^ ~ N/A for N divisible by 2, and epMe ~ iV/16 for N divisible by 4. For 
larger p, the mixing time tends to the classical value of iV^/16e. Although for 
N divisible by 4, the (coin-decohered) mixing time shows a minimum below the 
classical value at p ~ 2/N, this mixing time is still quadratic in iV. Thus although 
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Fig. 5. Numerical data for mixing times on cycles of size N = 29 and A'^ = 30 (□), for 
coin (dotted), particle (dashed) and both (solid) subject to decoherence, using s = 0.01. 
Also N = 28 ( X ) for coin. Both axes logarithmic 

noise on the coin causes the even- TV cycle to mix to the uniform distribution, it 
does not produce a significant speed up over the classical random walk. 

For decoherence on the particle position, withp < 16/7V^, epM^ ~ l/{N/2 — 
1) for TV divisible by 2, and spM^ ~ l/(-/V/4 + 3) for N divisible by 4. At 
p ^ Ib/N^, there is a minimum in the mixing time at a value roughly equal 
to the (A'' + l)-cycle pure quantum mixing time, mJ™'"-* ~ aN/s (with a a 
constant of order unity). Decoherence on the particle position thus causes the 
even-A^ cycle to mix in linear time for a suitable choice of decoherence rate 
p(min) ^ 16/A'^, independent of e so long as £ < 1/A'. 

For all types of decoherence, the odd-A^ cycle shows a minimum mixing time 
at a position somewhat earher than the even-A^ cycle, roughly p = but 
because of the oscillatory nature of P{x,p,T), the exact behaviour [25] is not a 
smooth function of p or e. As decoherence on the particle (or both) increases, 
at j3 ~ 16/ N'^, the mixing time passes through an inflexion and from then on 
behaves in a quantitatively similar manner to the adjacent-sized even- AT cycles, 
including scaling as Aii™™-* ~ aN/e at the inflexion. Thus for at least < p < 
the mixing time stays linear in A^, and exhibits the quantum speed up 
over the classical N'^. 

2.3 Decoherence in a Quantum Walk on a Hypercube 

Recalling from Sect. 1.6, we arc interested in the hitting time to the opposite 
corner, which was shown by Kempe [18] to be polynomial, an exponential speed 
up over a classical random walk. Kempe discusses two types of hitting times, 
one-shot, where a measurement is made after a pre-determined number of steps. 
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Fig. 6. Hitting probability on a 9-dimcnsional hypcrcubc for one-siiot (left) and con- 
current (right), perfect walk (circles) , with p = 0.05 (dotted), p = 0.1 ~ 1/9 (solid). 
Classical hitting probability barely visible (dashed) 

and concurrent, whore the desired location is monitored continuously to sec if 
the particle has arrived. In each case, the key parameter is the probability Ph of 
finding the particle at the chosen location. Here we consider only target locations 
exactly opposite the starting vertex. We have calculated Ph numerically following 
the scheme of [18, 26] with a iV-dimensional coin and the Grover coin operator, 
defined in (19). Figure 6 shows how Ph is affected by decoherence. All forms 
of decoherence have a similar effect on Ph, reducing the peaks and smoothing 
out the troughs. For the one-shot hitting time this is useful, raising Ph in the 
trough to well above the classical value, so it is no longer necessary to know 
exactly when to measure. For p < 1/N, the height of the first peak scales as 
Ph{p) = -Pft(O) exp{— (iV -|- a)p}, where < a < 2 depending on whether coin, 
particle or both are subject to decoherence. Thus Ph decreases exponentially 
in p, but p ~ 1/N only lowers Ph by a factor of 1/e, still exponentially better 
than classical. Continuous monitoring of the target location as in the concurrent 
hitting time is already a type of controlled decoherence, no new features are 
produced by the addition of unselectivc decoherence, but there is still a range of 
< p < 1/N within which the quantum speed up is preserved. 

2.4 Summary of Decoherence Effects in Coined Quantum Walks 

For the walk on a line, whilst the effect of decoherence is to reduce the standard 
deviation rapidly towards the classical value, it is possible to generate iiighly 
uniform distributions over a range proportional to T for small values of the de- 
coherence parameter, p. Uniform sampling is one of the basic tasks for which 
classical random walks are used, so being able to do this over a quadratically 
larger range with a quantum walk is certainly promising, though no quantum al- 
gorithms using this have been described to date. For a walk on a cycle the effects 
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of decoherence are more beneficial. An optimum rate of decoherence exists for 

which the rate of mixing is enhanced beyond the pure quantum bound. Further, 
any amount of decoherence removes the effect of the coin flip operator on the 
steady state, allowing all such walks to converge to a uniform distribution. Fast 
mixing to a uniform distribution is again an important basic property required 
for efficient random sampling, and is the limiting factor in many classical algo- 
rithms. For a hypercube, decoherence can still be tolerated so long its rate is 
kept smaller than 0{1/N), which is logarithmic in the system size (2^). 

Thus, for both experiments and algorithms, decoherence need only be con- 
trolled down to finite low levels, rather than negligible levels, in order to observe 
the intriguing quantum effects displayed by coined quantum walks. Many open 
problems remain concerning the best ways to do this, how to exploit the full 
power of the extra degrees of freedom provided by the quantum coin, and how 
to make quantum walks perform useful tasks over more general graphs that 
provide real quantum computational advantages over classical algorithms. 
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